Fluctuating flow angles and anisotropic flow measurements 
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Event-by-event fluctuations in the initial density distributions of the fireballs created in relativis- 
tic heavy-ion collisions lead to event-by-event fluctuations of the final anisotropic flow angles, and 
density inhomogeneities in the initial state cause these flow angles to vary with the transverse mo- 
mentum of the emitted particles. It is shown that these effects lead to characteristically different 
transverse momentum dependencies for anisotropic flow coefficients extracted from different exper- 
imental methods. These differences can be used to experimentally constrain flow angle fluctuations 
in the final state of heavy-ion collisions which, in turn, are sensitive to the initial state density 
fluctuations and the shear viscosity of the expanding fireball medium. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 



Due to quantum fluctuations of the positions of the 
nucleons inside the colliding nuclei, and of the positions 
of the colored quark and gluon constituents inside each 
nucleon, the density of the fireball matter created in colli- 
sions between ultra-relativistic heavy ions is highly inho- 
mogeneous in the plane transverse to the beam direction 
and fluctuates from event to event, even for collisions 
with identical impact parameters. Experiments at the 
Relativistic Heavy Ion Collider (RHIC) [1] and their the- 
oretical interpretation have established that this matter 
quickly thermalizes into a quark-gluon plasma (QGP) 
[2]. The initial density inhomogeneities lead to highly 
anisotropic pressure gradients, causing an anisotropic col- 
lective expansion of the fireball whose harmonic flow co- 
efficients v n and associated flow angles ty n (both defined 
below) fluctuate from collision to collision [3]. 

While v n fluctuations and the effect of their variance 
on different methods for measuring v n have been studied 
extensively over the last few years (for recent reviews and 
references to the original literature see [4, 5] ) , flow angle 
fluctuations and correlations have only recently found at- 
tention [6-14]. Gardim et al. [14] pointed out that, since 
the fluctuating flow angles \&„ depend on transverse mo- 
mentum pt and rapidity y, the usually assumed (and 
experimentally observed [15-18]) factorization of the az- 
imuthal oscillation amplitudes of the two-particle angular 
correlations into a product of single-particle flow coef- 
ficients is slightly broken even if these correlations are 
entirely due to collective flow. We show here that the 
PT-dependence and fluctuating nature of the flow angles 
ty n also affects the pr-dependence of the experimentally 
measured differential flow coefficients v n {jPT)i and that it 
does so in different ways for different experimental meth- 
ods of determining v n (j>T)- 
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Here <f> is the azimuthal angle around the beam direction 
of the particle's transverse momentum px, and the curly 
brackets denote the average over particles from a single 
collision. 1 Eq. (1) defines the flow coefficients and asso- 
ciated flow angles for the entire event, whereas Eq. (2) 
is the analogous definition for the subset of particles in 
the event with a given magnitude of the transverse mo- 
mentum px- We suppress the dependence of both types 
of flow coefficients on the rapidity y. v n are known as 
the "integrated" anisotropic flows, v n (px) are called "dif- 
ferential" flows. By definition, both v n and v h {pt) are 
positive definite. Hydrodynamic simulations show that 
in general the flow angles \& n depend on px, and that, as 
a function of px, ^(pt) wanders around the "average 
angle" ^ n that characterizes the integrated flow v n of the 
entire event (see Fig. 1 below and also Fig. 2 in Ref. [13]). 
Some theoretical and experimental definitions of vi have 
yielded values that turn negative over certain px ranges; 
we will see that this is due to defining the flows of each 
event relative to a fixed azimuthal angle (for example, 
relative to the direction of the impact parameter of the 
collision in theoretical calculations, or relative to the "in- 
tegrated" elliptic flow angle '5 2 m experiment), and that 
the same thing can happen for higher order harmonic 
flow coefficients when defining them relative to a fixed 
(i.e. px-independent) flow angle ^f n . The subject of this 



For each collision event the momentum distribution of 
finally emitted particles can be characterized by a set of 
harmonic flow coefficients v n and flow angles ty n through 



1 The average can include all charged particles or only particles of 
a specific identified species; we will not clutter our notation to 
account for these different possibilities. 
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FIG. 1: (Color online) (a) pt dependence of the flow angles ^ n (pr) relative to their average value 9 n (n= 1, . . . , 6), for pions 
from a single but typical Pb+Pb collision event at LHC energies in the 20—30% centrality class, computed with the viscous 
hydrodynamic code VISH2+1 using an initial density profile from the Monte-Carlo-Glauber model and r)/s = 0.08. Note the 
flip of the dipolar (n = 1) flow angle from to — n (happening around pt ~ 0.9 GeV in this particular event); this reflects the 
constraint from total transverse momentum conservation [19] . (b,c) The variance a of the cosine of the flow angle fluctuations 
(which is free of ambiguities of the angles by multiples of 2ir/n) as a function of pt, for pions (b) and protons (c) from 11,000 
hydrodynamic events in the same centrality class from which the event in (a) was taken. The variance is minimal around 
Pt — 1 GeV for pions and around pt = 1.5—2 GeV for protons (except for the dipolar flow angle which flips by 7r exactly in that 
momentum region, see panel (a)). It is smallest for the elliptic flow angle fluctuation >&2(pt) — $2, and becomes large for all 
flow harmonics at low pt- 



paper is to elucidate the origins of such differences be- 
tween different anisotropic flow measures and, in partic- 
ular, the manifestation of event- by- event fluctuations of 
the pT-dependent flow coefficients v n (pr) and flow angles 
V? n (pT) in different experimental flow measures. 



II. DIFFERENTIAL FLOWS FROM THE 
EVENT-PLANE METHOD AND FROM 
TWO-PARTICLE CORRELATIONS 

The key experimental difficulty is that, due to the fi- 
nite number of particles emitted in a each collision, the 
left hand sides of Eqs. (1,2) cannot be determined ac- 
curately for a single event. The V n are characterized by 
probability distributions that depend on the studied class 
of events (system size, collision energy and centrality) 
from which each collision takes a sample. Experimental 
flow measurements rely on a number of different methods 
that amount to taking different moments of that prob- 
ability distribution, by averaging over large numbers of 
events. Understanding the nature of these moments and 
reconstructing them from theoretical event- by- event dy- 
namical simulations are essential steps in a meaningful 
comparison between theory and experiment. 

Our main interest lies in the event- by- event fluctua- 
tions in the initial state of the collision fireball. These 
are primarily caused by the finite number of nucleons (or 
effective collision centers) in the colliding nuclei and un- 
related to detector capabilities. In addition, there are 
fluctuations related to the finite number of particles pro- 
duced (or detected) in the event which depend on colli- 
sion energy and (in part) on detector capabilities. They 
reflect the fact that in practice the final state of the fire- 



ball evolution, which in principle (with the appropriate 
dynamical evolution model) can be predicted from the 
initial state with perfect precision, cannot be measured 
with perfect precision, due to finite sampling statistics. 
In this paper we are not interested in the fluctuations 
arising from finite sampling statistics. The consequences 
of finite number statistical effects in a single event on the 
new observables proposed in this paper deserve a care- 
ful investigation that we leave for the future. Here we 
focus on the hydrodynamical consequences of unavoid- 
able event-by-event fluctuations in the initial state over 
which we have no control since they are rooted in the 
internal structure of the colliding nuclei, and with which 
we therefore have to live in any case even after we cor- 
rect the measurements for finite final state multiplicity 
effects. 

The most extensively used experimental methods for 
measuring anisotropic flows are the event-plane and two- 
particle correlation methods [4] . We begin with a discus- 
sion of the latter. Two-particle azimuthal correlations re- 
ceive contributions from the anisotropic collective flow as 
well as from non-flow correlations; the latter can be mini- 
mized by appropriate experimental cuts and corrected for 
[4, 20]. Again, we are not interested in non-flow correla- 
tions and will here simply ignore their existence, assum- 
ing that they have been corrected for in the experimental 
analysis. 

Two-particle correlation measures of anisotropic flow 
are based on correlators of the type 

({e'"(*- fc '}) (3) 

where <f>\ and <p2 are the azimuthal angles around the 
beam direction of two particles with transverse momenta 
Pti and px2, and (...) denotes the average over N cv ^> 1 
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events from a set of given characteristics (e.g. of collisions 
in a certain centrality bin), 



JV e - 



(0) = ({0}) := F a°}» 



(4) 



i=l 



whereas {...}, is the average of the observable O over all 
(or a specified subset of all) particle pairs in the event i: 
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Different chosen subsets for the event- wise average {•••}* 
define different correlation measures for the anisotropic 
flow coefficients as we will explain below. We note that 
throughout this paper we will always correlate pairs of 
particles of the same kind (e.g. protons with protons or 
charged hadrons with charged hadrons, but not protons 
with charged hadrons), unless specifically stated other- 
wise. We will also assume that they have the same rapid- 
ity y; generalization to particles with different rapidities 
is straightforward, following the procedure discussed be- 
low when we go from particles with the same to particles 
with different pr- 

The magnitudes v n (pT) of the anisotropic flow coeffi- 
cients defined in Eq. (2) fluctuate from event to event 
according to some probability distribution P(v n (pT))- 
Let us denot e the rm s mean of this distribution by 
t>„[2](pr) := \J (v^(pt)), and similarly the rms mean for 
the integrated flow v n by v n [2] :— (v 2 .). These rms 
means can be obtained from two-particle correlators of 
the type (3) as follows: 



«n[2](Pr) = ({e m ^-^} PT ) = ({e m ^} PT { 
t£[2] - ({e^** 1 -*■>}) = {{e in ^}{e- in * 2 }). 



PT/i 

(6) 



Note that for the differential flow in the first line of equa- 
tion (6), both particles are taken from the same pt bin, 
and that the event-wise pair averages {e 1 "^ 1 "^ 2 )} fac- 
torize in each event due to our assumptions (absence of 
non-flow two-particle correlations, independent hydrody- 
namic emission of particles 1 and 2). Due tolf>2 sym- 
metry under particle exchange, the exponential can be 
replaced by the cosine, and we get 



v 2 „[2](Pt) = 
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where A(f> = <f>\~(\)2 and the pair distribution has already 
been integrated over the average angle </>= (0i+</>2)/2. 

Note that in Eqs. (7) the single-event averages are nor- 
malized by the number of pairs in the event, before av- 
eraging over events. This is important: Since the pair 



multiplicity fluctuates from event to event and within 
a multiplicity bin, and multiplicity anti-correlates with 
impact parameter with which the magnitudes of some 
of the anisotropic flow coefficients are geometrically cor- 
related, this event-wise normalization avoids biasing the 
measured flow coefficients towards their values in events 
with larger than average multiplicity. 

Our definition of the integrated flow v n [2] agrees with 
the standard definition for the "two-particle cumulant" 
flow v n {2} [21-23], but the same is not true for the differ- 
ential flow t> n [2](pr) which differs from v n {2}(p^). The 
experimental definition of v n {2}(p^) is [21-23] 



v n {2}( PT ) := ({e^} PT1 {e-^})K{2} 

'vn(prKcos[n(tf n (pr)-*n)]V«n[ 2 ] . (8) 



Here only the first of the two particles within an event 
is taken from the desired px bin and particle species; 
it is correlated with all other particles detected in the 
event, with obvious statistical advantages compared with 
v„[2](pt) which requires both particles to be of the 
same kind and from the same px bin. The normaliza- 
tion factor is the total rms flow of all charged hadrons. 
The last expression shows that w„{2}(px) reduces to 
v n [2](p T ) = V ( w «(pt)) if and only if the flow angle 
does not depend on px, the event-by-event fluctuations 
of v„(px) affect only its normalization but not the shape 
of its pt dependence, and the v n fluctuations of the par- 
ticle species of interest are proportional to those of all 
hadrons. All of these assumptions are violated in hydro- 
dynamic simulations of bumpy expanding fireballs. The 
difference between v„{2}(px) and v„[2](px) is thus sen- 
sitive to event- by- event fluctuations of the px-dependent 
difference v P n (PT)~^'n between the flow angle of particles 
with momentum p^ and the average event flow angle, 2 in 
addition to the (largely independent) fluctuations in the 
magnitudes of v n and w„(pt)- 

Another approach to isolating effects arising from the 
px-dependence of the flow angles is a comparison of the 
PT-dependent rms flow v„ [2] (px) with the so-called event- 



2 More precisely, u„{2}(pt) depends on the difference between 
*n(PT) of the particles of interest and the average flow angle 
ty n of all detected particles. We checked numerically that the 
average hydrodynamic flow angles ^ n for identified pions and 
protons agree with great precision with the average flow angles 
for all particles in the event: Computing the ensemble average 
of (cos[n(*^— for all harmonics n and all collision central- 
ities, we found deviations of less than 1—2% in all cases except 
for some of the high-order harmonics with n > 6 whose calcula- 
tion is plagued by numerical errors at low px^0-2GeV arising 
from the finite grid spacing of our square numerical grid used in 
solving the hydrodynamic equations. 
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The equality in the second line arises from Eq. (2). Here 
for each event the "average flow angle" ^ n is first ob- 
tained by computing the Q n vector [24] 



Qn — Qn^ 



k=l 



(10) 



(where N is the number of detected particles in the event) 
and determining its phase. In principle, different choices 
for the weights uik can be considered [4], but for consis- 
tency with Eq. (9) one must choose Wfe = l. The "aver- 
age angle" ^ n for the event extracted from Q n in gen- 
eral depends on the types of particles included in the 
sum in Eq. (10). As noted in footnote 2, however, the 
average flow angle for particles emitted from a hydro- 
dynamic source is (within numerical precision) the same 
for all particle species, and the precision of extracting \I>„ 
in experiments can thus be maximized by including all 
detected particles in the Q n vector (10). 5 

The last line in Eq. (9) makes it clear that the differ- 
ential event-plane flows i>„{EP}(px) are sensitive to the 
event-by-event fluctuations of the px-dependent flow an- 
gles \&„(pt) around the "average flow" angle \P„. Just 
like the finite number statistical fluctuations 6 of the flow 



3 One can replace the cosine function in this definition by the expo- 
nential, omitting taking the real part in the second line, since the 
flow-angle fluctuations are symmetrically distributed such that 
the imaginary part vanishes after taking the event average (this 
has been verified numerically). 

4 Note that we define the n th -order event-plane flow relative to 
the ra th -order flow plane ty„, and not relative to the elliptic flow 
plane ^2 as sometimes done. 

5 Since in this paper we ignore finite particle statistical fluctuations 
in the final state, we know ty n with infinite precision for each 
particle species, and we will simply use these particle-specific 
values in our numerical results below. 

6 Due to the finite number of particles detected in each event, the 
accuracy of determining is limited by finite number statistics, 
and an accurate experimental estimation of the event-plane flow 
t> n {EP} requires an "event-plane resolution correction" [4]. As 
shown in [3, 20] (see also the discussion in [14]), which moment of 
the underlying v n distribution is actually measured by the total 
event-plane flow t)„{EP} depends on this event-plane resolution: 
for perfect resolution v n {EP} approaches the average flow (v„) 
whereas in the case of poor resolution it is closer to the rms 
flow v n [2] =v n {2} [3, 20]. The mathematical analysis in [20] 
applies only to the integrated flow which allowed to ignore the 
Px-dcpendcncc of v n fluctuations as well as initial-state related, 
PT-dcpcndcnt fluctuations of the flow angles that are not caused 
by finite multiplicity in the final state. In view of the latter, 
event-plane resolution effects on differential flow measurements 
and their correction require a new analysis. 



angle reconstructed from Q n around the "true" flow an- 
gle of the event, these fluctuations smear out the az- 
imuthal oscillations of the transverse momentum spectra 
and thus reduce the oscillation amplitudes w„{EP}(pr)- 
In contrast to the former, they arise from fluctuations in 
the initial state and thus cannot be eliminated by improv- 
ing or accounting for the resolution of the measurement 
of the final state. They carry valuable physical infor- 
mation about the initial state and the dynamics of its 
evolution into the final state. 

We can remove the sensitivity of the measured quantity 
to the j»T-dcpendcnt fluctuations of the flow angle by 
first computing for each event the magnitude v u (pt) of 
{e m *} PT = Vnipr), before summing over events: 



(Mpt)) = (\{e m % T e-^\) = {\{e m % T \) 



^/{cos(n0)}2 T + {sm{n4>)} 



(11) 



Since the quantity inside the event average does not de- 
pend on the average flow angle \I>„, this observable is not 
subject to an event-plane resolution correction. However, 
due to finite multiplicity in the final state, the right hand 
side will still in general be positive and non-zero experi- 
mentally even if there is no underlying anisotropic flow in 
the event. Again, how to properly account for such finite 
sampling statistical effects requires additional analysis. 

By comparing (v n (p T )) (11) with v n [2](p T ) (6,7), 
v„{2}(px) (8), and w„{EP}(px) (9), we can experimen- 
tally assess and separate the relative importance of event- 
by-event fluctuations in the magnitudes and directions of 
the anisotropic flows as functions of pr- 

Let us now proceed to two-particle correlations be- 
tween particles of different (but specified) momenta. 
Since in the first line of Eq. (7) both particles are taken 
from the same bin in pt, the flow angle & n (pr) drops 
out from the expression. This is not true for azimuthal 
correlations between two particles with different px [14]. 
In this case one finds [14, 25] 

Ka(pti,Pt 2 ) := ({e mWl " 02) } PTlPT2/ 
= <{e m01 We- m02 } PT2 ) = (V n (p Tl )V,: , . 
'v n (p T1 )v n (p T2 )e in ^^-*^») 



= [v n {pTi)v n (pT2)cos[n{^ n {p T1 )-^ n (p T2 ))}y (12) 

Due to parity symmetry, V n A(PTi,PT2) is real: while the 
quantity inside the event average (...) is in general com- 
plex for each individual event, its imaginary part averages 
to zero when summed over many events. 

To properly account for multiplicity fluctuations, in 
Eq. (12) the averages {. . . } PTi within an event are once 
again normalized by the total number of particles in- 
cluded in the average, similar to Eq. (7). For this 
reason, Kia(pti,PT2) = ({cos(nA</>)}p TliPT2 ) defined in 
Eq. (12) is not identical with the experimental quan- 
tity V n A(pTi,PT2) which is obtained from a Fourier de- 
composition with respect to the difference angle A(p 
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of the two-particle distribution obtained by summing 
over many events, without normalizing the contribution 
from each event by the corresponding event multiplic- 
ity [15-18, 25]. For a meaningful comparison between 
theory and experiment, one should either normalize on 
the experimental side the contribution from each event 
to the two-particle distribution by the number of pairs 
in the event, or weight the theoretical prediction for 
Wn(pTi)wn(pT2)e m( *" (pTl) "*" (PT2)) for each event i with 
a factor N^ irs / (N pa _i rs ) before summing over events. We 
prefer the first option since it avoids the geometric bias 
arising from the correlation between collision geometry 
and particle multiplicity. 

Equation (12) makes it obvious that the two-particle 
correlation coefficient V^a(pti,PT2) does not factorize 
into a product of single-particle anisotropic flow coeffi- 
cients [14]. There are two contributions to this break- 
ing of factorization: p-p-dependent event-by-event fluctu- 
ations of the magnitude of the flow coefficient v n , and 
px-dependence of the flow angles \I/„ [14] (which also 
fluctuate from event to event). It is possible to define 
a non-factorizing correlator that is only affected by the 
fluctuations of Vu(pt) but insensitive to the flow angles: 

= ( v /{cos(nA0)}2 Ti pT2 + {sin(nA0)}2 Ti pT2 ) . (13) 

It is obtained experimentally by first obtaining the mag- 
nitude of the quantity {e m ^ 1 } PT1 {e _m * 2 } J , T2 for each 
event, normalizing it to the number of pairs used for 
its computation, and than adding the results for many 
events. Its sensitivity to finite number statistical ef- 
fects should be similar to Eq. (11) and needs to be ex- 
plored. By comparing the quantity V^a(pt1)Pt2) from 
Eq. (12) with (v n (pTi)v n (pT2)} from Eq. (13) one can 
assess the importance of the pr-dependence and event- 
by-event fluctuations of the flow angles ^ n (which affect 
the former but not the latter). 

III. THE EFFECT OF FLOW FLUCTUATIONS 
ON DIFFERENTIAL v n MEASURES 



the final stable hadron spectra. The details of the hy- 
drodynamic simulations, initial conditions and freeze-out 
parameters are not important for the qualitative study 
presented here, but the interested reader can find them 
described in Refs. [12, 26, 27]. Here we only note that 
MC-Glauber (MC-KLN) initial conditions were hydrody- 
namically evolved with specific shear viscosity rj/s = 0.08 
(0.2). 

We present results for pions and protons, represent- 
ing light and heavy particle species. Qualitatively, al- 
though not quantitatively, the same generic features are 
observed with MC-KLN and MC-Glauber model initial 
density profiles, and we show examples of both. Fig- 
ures 2 and 3 show elliptic and triangular flows in their 
left and right panels, for central (Fig. 2) and peripheral 
(Fig. 3) Pb+Pb collisions. The reader should compare 
the curves for t> 2 .3[2] and (1*2,3)5 which are not affected 
by flow angle fluctuations (c.f. Eqs. (6) and (11)), with 
those for u 2> 3{2} and w 2 ,3{EP}, which are affected by the 
px-dependence of the flow angles \Jj„ and their event-by- 
event fluctuations (c.f. Eqs. (8) and (9)): For protons 
with transverse momenta below about 1 GeV, flow angle 
fluctuations are seen to cause a significant suppression of 
the latter (in some cases even leading to negative elliptic 
flow values). 7 For the much lighter pions flow angle fluc- 
tuation effects are almost invisible at low • For protons 
they gradually disappear, too, as one goes from central 
(Fig. 2) to peripheral (Fig. 3) collisions. 8 

Event-by-event fluctuations of the magnitudes of 1*2,3 
are accessible by comparing (1*2,3) with 1*2,3 [2] = (v% s) 1 ^ 2 - 
When plotting the ratios V2,3[2](pt)/(«2,3)(pt) for cen- 
tral (0—5% centrality) collisions, where anisotropic flows 
are caused exclusively by fluctuations, with negligible 
geometric bias from a non-zero average deformation of 
the nuclear overlap region, we found for both pions 
and protons a constant (i.e. px-independent) value of 
2/Vtt«1.13. This is expected [28, 29]: If the flow an- 
gle \I/„ is randomly distributed relative to the reaction 
plane, the components of V n (pr) along and perpendicu- 
lar to the reaction plane are approximately Gaussian dis- 
tributed around zero, and the magnitude u„(pt) of the 
complex flow coefficient is Bessel-Gaussian distributed 



In this section we compare the differential flows v u (pt) 
extracted from the 22,000 viscous hydrodynamic simula- 
tions per centrality bin of 2.76 A TeV Pb+Pb collisions at 
the LHC (11,000 each with MC-Glauber and MC-KLN 
initial density profiles) that were generated in Ref. [12]. 
We use the Cooper-Frye prescription to compute from 
the hydrodynamic output on the freeze-out surface the 
single-particle distributions dN '/ '(dypTdprdcj)) as contin- 
uous functions of pr and <p (i- e - wc do not sample the 
distribution to generate a finite number of particles per 
event, but pretend that the spectrum is sampled infinitely 
finely - this avoids the need to correct for effects arising 
from finite number statistics, such as imperfect event- 
plane resolution). All resonance decays are included in 



7 Note that the factor cos[n(^ n {pr) — *„)] in Eqs. (8) and (9) is 
maximal if 'I'n(PT) * s always aligned with The suppression 
of, say, v„{EP}(pt) relative to (v n )(pT) does therefore not indi- 
cate a definite momentum tilt of the emitting source at a given 
Pt relative to the average but simply reflects a nonzero 
difference ^ n (pT)—^n that fluctuates from event to event, sup- 
pressing the value of cos[n( 1 I'„(px) — *n)] for cither sign of the 
difference. 

8 The curves shown in Figs. 2, 3 include the decay products from 
unstable hadronic resonances. We have observed that for protons 
the flow angle fluctuation induced difference at low pp between 
( v 2, 3[2], (^2,3)) on the one hand and (i>2,3{2}, «2,3{EP}) on the 
other hand doubles if only directly emitted ( "thermal" ) particles 
are included in the analysis. Resonance decays thus dilute the 
sensitivity of the proposed observables to flow angle fluctuations 
by about 50%. 



6 




FIG. 2: (Color online) Comparison between the different definitions of the differential flows d„[2](pt) (6,7), « n {2}(px) (8), 
v n {EP}(pT) (9), and («n(pT)) (11), for pions and protons from central (0—5% centrality) Pb+Pb collisions at *Js = 2.76 A TeV, 
computed with the viscous hydrodynamic code VISH2+1. See text for discussion. 




FIG. 3: (Color online) Same as Fig. 2, but for peripheral Pb+Pb collisions at 60—70% centrality. 



with V( v Upt)) = ^(UnOr)) ( see Eqs. (4) and (5) in 
Ref. [28]). A similar px- independent ratio is not ob- 
served at larger impact parameters: even for triangular 
flow, which continues to be fluctuation-dominated also 
at non-zero impact parameters, we observe deviations of 
the ratio i>3[2](pt)/(^3)(pt) from 2/y/n at both low and 
high px; for elliptic flow these deviations are larger and 
significant at all px. 

Interestingly, for central collisions we found approx- 
imately the same constant value 2/y / 7r for the ra- 
tio W2,3{2}(pt)/w2,3{EP}(pt) (except near the px val- 
ues where either the numerator or denominator passes 
through zero). Looking at the definitions (8) and (9), 
this suggests an approximate factorization of the px- 
dependent flow angle fluctuations (which enter through 
the factor cos[n(^ n (pT)—^ n )) that cancels between nu- 
merator and denominator if it fluctuates independently) 
from the fluctuations of the magnitude v„(px), as well as 
an approximate px-independence of the V2.3 fluctuations. 

To follow up on these observations and gain deeper 



insight into the relative importance of flow angle fluctu- 
ations in different px ranges, let us look at Figs. 2, 3 and 
note that the frequently measured quantity t>2,3{2}(px) 
behaves like the event-plane flow ^2,3{EP}(px) at low 
Px and like the differential rms flow ^2,3 [2] (px) at in- 
termediate px- This suggests that it is dominated by 
flow angle fluctuations at low px and by fluctuations 
of the magnitude of i>2,3(pr) at higher px). In central 
collisions, the proton W2,3{2}(px) even turns negative at 
low pt, whereas 1*2,3 [2] (px) is by definition always pos- 
itive. A related observation is that the proton event- 
plane flow i>2 3 {EP}(pt) m Eig. 2 approximately agrees 
with v| 3 {2}(pt) at low px (where flow angle fluctua- 
tions seem to have strong effects) but with the mean flow 
(^2 3 (pt)) at higher px (where flow angle fluctuation ef- 
fects are weak). This is reminiscent of the behavior of 
the pt -integrated event-plane flow which approaches the 
mean flow for good event-plane resolution. Flow angle 
fluctuations appear to have similar effects on flow mea- 
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FIG. 4: (Color online) (a,b) Ratios of differently defined elliptic (a) and triangular (b) flow coefficients for pions and protons as 
functions of pt- (c,d) px-dependence of the separate fluctuations of the magnitudes v n and the angles ty n for pions and protons 
as discussed in the text. All curves from viscous hydrodynamics with ?7/s = 0.08 for central 2.76 A TeV Pb+Pb collisions with 
MC-Glauber initial conditions. Results for MC-KLN initial conditions evolved with n/s = 0.2 look very similar. 



sures as a decrease in flow angle resolution. The differ- 
ence is that the former is a physical effect due to initial- 
state fluctuations whereas the latter is a finite sampling 
statistical effect in the final state and affected by detector 
performance. 

To make these qualitative observations quantitative, 
we plot in the upper two panels of Fig. 4 the ratios 
U2,3{2}/u2,3[2] and U2,3{EP}/(«2,3) as functions of pr, for 
both pions and protons. (We focus here on the results 
from Fig. 2 for central collisions where all anisotropic 
flows are fluctuation-dominated.) In each case the nu- 
merator is sensitive to the flow angle fluctuations while 
the denominator is not. However, numerator and de- 
nominator are also differently affected by fluctuations 
in the magnitudes of v n . Both ratios are seen to be- 
have very similarly, staying close to 1 at intermediate 
Pt but dropping steeply at low pt and more moderately 
at high pt- The steep drop at low transverse momenta 
sets in at pt ~ 1 GeV for protons, but at much smaller 
Pt < 0.25 GeV for pions. We do not have a full under- 
standing of this mass dependence, beyond the qualita- 
tive observation that the minimum of the variance of the 
flow angle fluctuations shown in Figs. lb,c is shifted to 
higher pt for protons compared to pions, and that quite 
generally strong radial flow shifts all flow anisotropies to 
higher pt values for heavier particles. 

The lower two panels of Fig. 4 demonstrate that the 
behavior of the ratios shown in the two upper panels 



is strongly dominated by flow angle fluctuations. The 
dashed lines in Figs. 4c, d show the flow angle fluctuations 
(cos[n(^f n (pT)—^ n )]) in isolation. Their pt dependence 
alone is almost sufficient to completely explain the shape 
of the curves in panels (a) and (b). The solid lines in 
Figs. 4c, d show that at intermediate pt fluctuations in 
the magnitudes of the pT-dependent flow v n (pr) and the 
px-integrated v n tend to be correlated with each other 
(v n (jpr) oc v n ) while they appear to fluctuate more inde- 
pendently at low and high pt- At high pt this decorre- 
lation contributes to the suppression of the ratios shown 
in panels (a,b). At low pr, the decorrelation of the px- 
dependent flow magnitude fluctuations v u (pt) from the 
px-integrated flow v n does not become effective until af- 
ter the ratios have already been suppressed by flow angle 
fluctuations, and its effect is therefore subdominant. 

In summary, we see for central collisions that at low 
Pt the differences between v„{2}(pt) and v„[2](pt), as 
well as between v„{EP}(pt) and (v n (pr)}, are dominated 
by flow angle fluctuations, whereas at high px fluctua- 
tions of both the flow angles and flow magnitudes must 
be considered to explain their differences. At intermedi- 
ate px flow angle fluctuations appear to be unimportant, 
v„(pt) fluctuates in sync with the pT-integrated v n , and 
the differences between i>„{2}(pt) and v„[2](pt), as well 
as between v„{EP}(pt) and (v n (pr)), vanish. 

Figure 5 shows the same ratios as Fig. 4 for peripheral 
Pb+Pb collisions, again using MC-Glauber initial con- 
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FIG. 5: (Color online) Same as Fig. 4, but for peripheral collisions at 60—70% centrality. 




FIG. 6: (Color online) Similar to Figs. 2,3, but with MC-Glauber initial conditions for Pb+Pb collisions in the 0—10% centrality 
bin. For clarity only f n {EP}(px) (Eq. (9), thick dashed) and (v„(pt)) (Eq. (11), thin solid) are shown, but for all harmonics 
from n = 2 to n — 9 (scaled by appropriate factors for best visibility). This set of plots focusses on the low-px region px < 1.5 GeV 
where the effects from flow angle fluctuations are strongest. See text for discussion. 



ditions with ij/s = 0.08. 9 Compared to central collisions 



9 The main difference with results from MC-KLN initial con- 
ditions with r]/s = 0.2 (not shown) is that the latter exhibit 
stronger suppression effects from the flow fluctuation factor 
cos[n('J' n (pT)— ^n)] in the high-px region px ^ 1 GeV (see also 
Fig. 3). 



(shown in Fig. 4), the flow angle fluctuation effects at 
low pt are much weaker and appear to be shifted to lower 
transverse momenta, for both pions and protons. At high 
Pt ^ 1 GeV, Figs. 5c, d show that effects from fluctuations 
of the flow angles (dashed lines) dominate over those from 
fluctuations of the flow magnitudes (solid lines). 

Finally, in Figure 6 we explore (for near-central colli- 
sions) how the flow angle fluctuation effects, which push 
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Pti (GeV) pti (GeV) p T1 (GeV) 

FIG. 7: (Color online) The ratios r2,3,4(pTi,PT2) (thick dashed lines) and f2,3,4(pTi,PT2) (thin solid lines), defined in Eqs. (14) 
and (15), as functions of pti for different pT2 ranges, as indicated. Filled symbols represent ALICE data for Pb-Pb collisions at 
\fs — 2.76 A TeV [16]. The lines are the corresponding viscous hydrodynamic calculations with MC-Glauber initial conditions, 
using specific shear viscosity r\js = 0.08. Panels (a,b,c) show r2,3,4 and r2,3,4 for central (0—10%) collisions, panels (d,e,f) show 
the same ratios for peripheral (40—50%) collisions. 



the event-plane flow w„{EP} at low-px below the value 
of the average flow (v n ), evolve as the harmonic order 
n increases. (For n > 4 we do not show results below 
Px = 0.2GeV, for technical reasons explained in foot- 
note 2.) For pions, flow angle fluctuations are invisi- 
ble in the shown px region for all flow harmonics; for 
protons, they are clearly visible for all harmonic flows. 
The relative magnitude of their effect on the difference 
(iVi)(pt)— v n {EP}(pT) at any fixed px decreases as n in- 
creases, but the difference remains nonzero over a larger 
Pt range for the higher harmonics. 



IV. NON-FACTORIZATION OF 
FLOW-INDUCED TWO-PARTICLE 
CORRELATIONS 

The breaking of factorization of flow-induced two- 
particle correlations by flow fluctuations was first em- 
phasized by Gardim et al. [14]. Their study was based 
on simulations using ideal fluid dynamics, which are here 
repeated with viscous fluid dynamics. A comparison of 
Figs. 7, 8 below with the plots shown in Rcf. [14] shows 
that viscous effects reduce the amount by which event-by- 
event fluctuations break factorization. We here explore 
the relative role played in this context by fluctuations in 
the magnitudes and angles of the flows. 



To this end we define the following two ratios, both 
symmetric in pxi an d PT2 '■ 



r n (PTl,PT2) 



Ka(pti,PT2) 



V n A (pTl , PTI ) V n A (PT2 , PT2 ) 
(V n ( P Tl)V*( PT2 )) 



V(\Vn( P TlW)(\V n (pT2W) 

(v n (pti ) v n (p T 2 ) cos [n(* n (pri ) - *n (PT2 ) ) ] } 



(14) 



V( V n(PTl))(v n (pT2)} 
f„(p T l,PT2) (15) 
,_ («n(PTl)^«(PT2)cOs[n(^»(p T1 )-^»(p T2 ))]) 
(v n {PTl)v n {PT2)} 

The ratio r n , first introduced and studied with ideal fluid 
dynamics in [14], is sensitive to fluctuations of both the 
magnitudes v n (pr) and angles $„(pt) of the complex 
anisotropic flow coefficients V n (pr) defined in Eq. 2). The 
second ratio f n , on the other hand, differs from unity only 
on account of flow angle fluctuations. By comparing the 
two ratios with each other and with experimental data we 
can isolate the role played by flow angle fluctuations in 
the breaking of factorization of the event-averaged two- 
particle cross section. In the absence of non-flow corre- 
lations both ratios are always < 1. 

Figures 7 show these ratios for all charged hadrons 
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FIG. 8: (Color online) Same as Fig. 6, but for MC-KLN initial conditions that have been evolved hydrodynamically with 
rj/s — 0.2 (a shear viscosity value that is 2.5 times larger than the one used in Fig. 6). 



as functions of pri <PT2 for fixed ranges of pr2, indi- 
cated by different colors. 10 Figs. 7a,b,c focus on cen- 
tral, Figs. 7d,e,f on peripheral collisions; in both cases, 
we used MC-Glauber initial conditions and evolved them 
with VISH2+1 using rj/s = 0.08 for the specific shear vis- 
cosity. In central collisions the hydrodynamic simulations 
appear to overpredict the factorization breaking effects, 
while in peripheral collisions theory and data agree some- 
what better. More precise experimental data would be 
desirable. The comparison of r n (dashed lines) with f n 
shows that a significant fraction ( ~ 50% or more) of the 
effects that cause the breaking of factorization arises from 
flow angle fluctuations. This seems to hold at all the 
transverse momenta shown in the figures. A comparison 
of the top and bottom rows of panels in Fig. 7 shows that 
factorization-breaking effects are stronger for harmonics 
that are fluctuation dominated (i.e. all harmonics in cen- 
tral collisions, and the odd harmonics (especially V3) in 
peripheral collisions) and appear to weaken for V2 and W4 
in peripheral collisions where both the magnitudes ^2,4 
and the flow angles ^2,4 are mostly controlled by collision 
geometry. 

To explore the effects of shear viscosity of the ex- 
panding fluid on the breaking of factorization we show 



The pT2 ranges are adjusted to the experimental data, and the 
ratios were computed by first averaging the numerator and de- 
nominator over the given px2 range. 



in Figure 8 the same data as in Fig. 7, but compared 
with hydrodynamic calculations that use MC-KLN ini- 
tial conditions evolved with ij/s — 0.2 (a 2.5 times larger 
viscosity than used in Fig. 7). Obviously, the MC-KLN 
model produces a different initial fluctuation spectrum 
than the MC-Glauber model, so not all of the differences 
between Figs. 7 and 8 can be attributed to the larger 
viscosity. However, in conjunction with the ideal fluid 
results reported in [14], the comparison of these two fig- 
ures strengthens the conclusion that increased shear vis- 
cosity tends to weaken the fluctuation effects that cause 
the event-averaged two-particle cross section to no longer 
factorize. 



V. SUMMARY 

All experimental precision measures of anisotropic flow 
in relativistic heavy-ion collisions are based on observ- 
ables that average over many collision events. It has been 
known for a while that both the magnitudes v n and flow 
angles ^ n of the complex anisotropic flow coefficients V n 
fluctuate from event to event, but only very recently it 
became clear that not only the v n , but also their asso- 
ciated angles ^ n depend on pr, and that the difference 
^nipr)— between the px-dependent and px-averaged 
flow angles also fluctuates from event to event. In the 
present study we have pointed out that these flow an- 
gle fluctuations leave measurable traces in experimen- 
tal observables from which the ensemble-averaged px- 
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dependent anisotropic flows are extracted. We have in- 
troduced several new flow measures and shown how their 
comparison with each other and with flow measures that 
are already in wide use allows to separately assess the im- 
portance of event-by-event fluctuations of the magnitudes 
and angles of V n = v n e m<S/n on experimentally determined 
flow coefficients. 

Viscous hydrodynamic simulations show that flow 
angle fluctuations affect the pr-dependent flow coeffi- 
cients of heavy hadrons (such as protons) more visi- 
bly than those of light hadrons (pions). In near-central 
collisions, where anisotropic flow is dominated by ini- 
tial density fluctuations rather than overlap geometry, 
the effects from flow angle fluctuations appear to be 
strongest for particles with transverse momenta pt^tu. 
A precise measurement and comparison of {v n (p^)) (Eq. 
(11)), v n {EP}(p T ) (Eq. (9)), v n [2]{p T ) (Eqs. (6,7)), and 
v„{2}(px) (Eq. (8) for identified pions, kaons and pro- 
tons with transverse momenta pr < 2 GeV should be per- 
formed to confirm the hydrodynamically predicted ef- 
fects from flow angle fluctuations. The theoretical in- 
terpretation of these measurements requires a reanalysis 
of finite sampling statistical effects on the p^ -dependent 
differential flows, stemming from the finite multiplicity 
of particles of interest in a single event, which we did 
not consider here. The proposed comparison holds the 
promise of yielding valuable experimental information to 
help constrain the distribution of initial density fluctu- 



ations in rclativistic heavy ion collisions and may prove 
crucial for a precision determination of the QGP shear 
viscosity. 

We also showed that flow angle fluctuations are respon- 
sible for more than half of the hydrodynamically pre- 
dicted factorization breaking effects studied in Ref. [14] 
and in Sec. IV above, and that these effects are directly 
sensitive to the shear viscosity of the expanding fluid, 
decreasing with increasing viscosity. By combining the 
study of various types of differential anisotropic flow mea- 
sures with an investigation of the flow-induced breaking 
of the factorization of two-particle observables into prod- 
ucts of single-particle observables one can hope to inde- 
pendently constrain the fluid's transport coefficients and 
the initial-state fluctuation spectrum. 
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